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ABSTRACT 

Two sequences of solar images obtained by the Transition Region and Coronal Explorer in three UV 
passbands are studied using wavelet and Fourier analysis and compared to the photosphcric magnetic 
flux measured by the Michelson Doppler Interferometer on the Solar Heliospheric Observatory to 
study wave behaviour in differing magnetic environments. Wavelet periods show deviations from the 
theoretical cutoff value and are interpreted in terms of inclined fields. The variation of wave speeds 
indicates that a transition from dominant fast-magnetoacoustic waves to slow modes is observed when 
moving from network into plage and umbrae. This implies preferential transmission of slow modes into 
the upper atmosphere, where they may lead to heating or be detected in coronal loops and plumes. 
Subject headings: Sun: chromosphere - Sun: magnetic fields - Sun: oscillations - Sun: UV radiation 



1. INTRODUCTION 

It has been believed for over forty years that oscilla- 
tory motion s play an important r ole in the outer solar 
atmosphere (jLeighton et al.Ml962|) . This is particularly 
true of the wide variety of perturbations that are possi- 
ble within the chromosphere, as this is the atmospheric 
layer which spans the transition from the domination of 
gas pressure to that of the magnetic field. These oscil- 
lations may be the signature of wave energy carried out- 
ward from the photosphere and deposited in the chro- 
mosphere, hence resulting in heating of the outer at- 
mosphere. As such, spatial variations in oscillatory be- 
haviour could yield information on the local topolog y 
of the magnetic env i ronmen t (e.g.. lFinsterle et al.l l20041. 
iFossum fc Carlssonl l)2005afl showed that hydrodynamic, 
high-frequency waves in the range 5 — 50 mHz do not 
contribute sufficient energy to heat the atmosphere in 
a regi on of weak magnetic field. Also, Soc as-Navarrol 
(2005) found that electric currents in the chromosphere 
above a sunspot do not show a substantial spatial cor- 
relation with regions of cither high temperature or large 
temperature gradient. Hence, it appears entirely plausi- 
ble that differing forms of magnetoacoustic waves could 
supply some portion of the excess energy input through 
their existence in both extremes of the chromospheric 
environment. 

Previous work on chromospheric oscillatory behaviour 
has reported on differences between spatial regions such 
as "cell boundary" and "cell interior" (e.g.. iDame et alJ 
l!984HDeubner fcFlecU 1990ft . or the (equally as ambigu- 
ously defined) magne tic "network" an d non-magnetic 
"inte rnetwork" (e.g., iLites et all 119931: iMcAteer et"aD 
|2004|) . Although the findings of these authors are gener- 
ally discussed in relation to the magnetic environment, 
spatial distinctions were often not based upon any mag- 
netic field measurements, but instead on the magnitude 
of time-averaged emission within some spectral line or 
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passband. This arbitrary form of classification may give 
the false impression that some sharp transition occurs 
when moving from the weak-field regime to that of higher 
field strengths. In reality, the observed distribution of 
solar magnetism shows a continuous variation of field 
strengths (Domingucz Cerdena et al.l l2006j) . suggesting 
that forms of oscillatory behaviour which are (in any 
manner ) magnetically related should vary continuously 
as well ({Lawrence et al.ll2003ft . 

The comparison of chromospheric oscillations to sur- 
faces of constant plasma-/?, as determined by potential- 
field extrapol ations from photospher i c longitudinal mag- 
netic fields llMcIntosh et all 120031: iMcIntosh fc~S milie 
12001 iFmsFerle et al.ll200|T has shown the spatial distri- 
bution of oscillatory power and propagation characteris- 
tics. We note that (3, the ratio of gas to magnetic pres- 
sure, is given by (2/ r y)(cs/vA) 2 , where cs is the sound 
speed (=7P/p), 7 is the ratio of specific heats (=5/3), 
P is the gas pressure, p is the mass density, va is the 
Alfven speed (=B / \JA-Kp), and B is the magnetic field 
strength. However, the above studies mainly report upon 
the co-spatiality of oscillation behaviour with a chosen 
plasma-/3 level (usually the unity surface; cs ~ va) with- 
out yielding qualitative diagnostic feedback on the local 
atmospheric environment. In addition, the complicated 
two-dimensional characteristics of chromospheric wave 
propagation, both hv drodvnamic ( [ Rosenthal et al.ll2002^l 
and magnetoacoustic l)Bogdan et al.ll2003|) in nature, are 
only now beginning to be successfully explored through 
careful modeling. The results of such work are decidedly 
sobering, indicating a veritable plethora of interacting 
wave modes propagating at different speeds, sometimes 
in differing directions. Another consideration apparent 
from the simulations is the importance of knowing the 
plasma environment (i.e., gas or magnetic domination) 
in which the waves exist. 

Our paper seeks to address the variable distribution 
of chromospheric UV intensity oscillations by studying 
their occurrence as a function of both detection height 
and the underlying photospheric magnetic flux. Obser- 
vational data and their sources are discussed in §[2] with 
descriptions of the reduction and alignment processes, 
while § outlines the wavelet and Fourier analysis. The 
time-localized nature of the wavelet transform is used to 



2 



D. S. Bloomfield et al. 



TABLE 1 

TRACE UV and mdi Data Set Summaries 



Data Set Date TRACE MDI 

Start Time Frames Cadence Initial Pointing Start Time Frames Cadence Initial Pointing 
(UT) W (UT) (s) 

1 2004 Jul 03 15:50 128 16.0 (6", -16") 15:37 1 (-10", -122") 

2 2000 Sep 22 08:00 312 22.7 (-174", 124") 08:01 118 60.0 (-88", 44") 



record the available distributions of wavepacket periods 
and durations. These parameters are of interest as oscil- 
lation periods can indicate which form of wave may exist 
(i.e., acoustic or magnetic) while the duration in terms of 
oscillatory cycles indicates the degree of energy carried 
if an oscillation is indeed the signature of a propagating 
wave. A Fourier phase difference technique is then used 
to directly study wave propagation characteristics. The 
results are presented and discussed in § 0] in terms of 
wave behaviours and the plasma environment, while in 
§ El the implications of our work are summarized. 

2. DATA REDUCTION 

The chromospheric data studied were obtained by 
the Transition Region and Coronal Explorer (TRACE; 
lHandv et aflll999]) . Two 256 x 256 arcsec 2 regions were 
recorded in the TRACE UV passbands with a sampling 
of 0.5 arcsec pixel -1 on 2004 July 3 and 2000 September 
22, with no compensation for solar rotation. Details of 
both data sets are given in Table^wbile Figure^depicts 
the time-averaged 1700 A emissions. In set 1, images 
were obtained in the order 1600 A, 1700 A, and 1550 A 
with a strict 4 s timing offset because of the TRACE fil- 
terwheel set up, resulting in an overall stable cadence of 
16 s for each individual passband, while set 2 comprised 
of successive 1550 A, 1600 A, 1700 A, and white-light 
images with an average cadence of 22.7 s. 

Initially, data from each UV passband were dark sub- 
tracted and flat fielded. Portions of any images with cor- 
rupted data were replaced by the temporal mean of the 
preceding and following images, and JPEG comp ression 
effects were removed. Following iDeForesn l)2004l) , data 
were derippled in both space and time by thresholding 
and removing spikes in the spatio-temporal Fourier do- 
main. A final sweep of despiking was applied to remove 
cosmic ray hits and persistent "hot pixels" were removed 
from the data. After this processing, each UV passband 
datacube was co-aligned by spatially cross correlating 
every image to the mid-image of the respective time se- 
quence. Although difficult to attribute heights of forma- 
tion (HOF) to the broadband UV emissions (due to the 
large range of temperature coverage involved in their pro- 
duction, and the large HOF overlap between the 1700 A 
and 1600 A passbands), we assumed that the HOF in- 
creased from the 1700 A passband to 1600 A and then t o 
1550 A ^McAteer et aU2004HFossurn fc Ca,rlssonl2005rl. 
Mid-images of the 1600 A and 1550 A data were sepa- 
rately co-aligned to the mid-image of the 1700 A dat- 
acube for each data set, to achieve alignment between 
passbands. The 1700 A data were taken as the refer- 
ence for the alignment process because this passband has 
the response function with the lowest formation height 
in the solar atmosphere out of the TRAC E passbands 
(form ed around the temperature minimum; Uudge et al.l 
|2001[K thus allowing for a better reference comparison for 
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Fig. 1. — Time— averae ed TRACE 1700 A images for set 
1 (upper) and set 2 (lower). Overlaid thin (thick) contours 
mark the mdi |50| G (|500| G) level. 



the alignment to the underlying magnetic field data. 

Photospheric line-of-sight (LOS) magnetic field data 
were obt ained by the M ichelson Doppler Interferome- 
ter (mdi; Sch errer et alJl!995[) on board the So lar He- 
liospheric Observatory (SoHO: fFleck et aTlll995|) . These 
data consist of longitudinal magnetic field measurements 
recorded with a spatial sampling of 0.6 arcsec pixel -1 , 
as MDI was operating in high-resolution mode. The MDI 
magnetograms were processed to Level-1.8 by the stan- 
dard data pipeline and compensated for the roll angle 
of the SoHO spacecraft as well as the difference in image 
scales due to the differing observing positions of the satel- 
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Fig. 2. — Upper: Frequency distribution of the total number 
of pixels (ordinate) in the selected mdi absolute field strength 
intervals (abscissa) for set f . Lower: Percentages of the total 
number of pixels which show detections in each of the three 
TRACE UV passbands. The distributions for detections in 
the 1700 A, 1600 A, and 1550 A data are shown in dotted, 
dashed, and dotted-dashed lines, respectively. 

lites. In the case of set 2, magnetograms were co-aligned 
in the same manner as the UV datacubes before taking 
the temporal mean. Finally, mdi pixels corresponding 
to the time-rotated TRACE field-of-view (FOV) were 
re-sampled to the TRACE pixel size. The correlation of 
the co-spatial portions of the two data sets is shown by 
the contours in Figure Q] This form of MDi-to- TRACE 
alignment was verified by spatially cross correlating the 
magnetic field magnitude to the intensity of the TRACE 
1700 A mid-image, achieving sub-pixel accuracy. 

It is important that the field strengths presented 
here are, as for any magnetogram instrument, inter- 
preted as the integra tion of net flux over a pixel area. 
iBerger fc Litesl l)2003fl have shown, through direct com- 
parison to the more reliable spectropolarimetry of the 
Advanced Stokes Polarimeter (asp), that mdi system- 
atically underestimates the true longitudinal flux when 
operated in full-disk mode. However, we choose to re- 
port raw MDI values because our work makes use of mdi 
data obtained in high-resolution mode and the Berger 
& Lites correction factors may not be accurate — the 
difference in pixel sampling, and thus spatial flux aver- 
aging, should affect the slope of the MDI-ASP flux-flux 
relation. 

A magnetic field binning scheme was used in order to 
study oscillatory parameters in the TRACE UV data as 
a function of the absolute MDI magnetic field. Initially 
intervals of equal size were used, but the number of pix- 
els decreases rapidly with increasing field strength. To 
counteract this, pixels with field strength values above 
1 20 1 G were instead separated into bins of systematically 
increasing size (listed in the first column of TableEJl ■ The 
lower limit to the mdi magnetic field strength intervals is 
directly based on the instrumental no ise level being re- 
ported as ~ 1 20 1 G ijScherrer et al.ll995lh The probability 
distribution function (PDF) of pixels used in our work 
is illustrated for the case of set 1 in the upper panel of 



Figure |2 while the lower panel depicts the percentage of 
pixels that contain wavelet oscillation detections. 

3. ANALYSIS METHODS 
3.1. Wavelet Power Transform 

One of the most prevalent tools that has been recently 
used in observational studies of oscillatory behaviour is 
the technique of wavelet analysis (see the excellent in- 
troduction provided by Torrence & Compo 1998 and its 
application to solar data in, e.g., De Moortel et al. 2000, 
McAteer et al. 2004, and Mcintosh & Smilie 2004). This 
multiscale analysis is better suited than Fourier analy- 
sis to both the spatial and temporal non-stationarity 
exhibited by oscillations in the solar atmosphere. For 
our study of oscillatory signals we chose the standard 
Morlet wavelet, a Gaussian-modulated sine wave which 
has an associated Fourier period, P, corresponding to 
1.03s, where s is the wavelet scale. The wavelet trans- 
form is a convolution of the time series with the analyzing 
wavelet function whereby the complete wavelet transform 
is achieved by varying the wavelet scale, which controls 
both the period and temporal extent of the function, and 
scanning this through the time series. Regions exist at 
the beginning and end of the wavelet transform where 
spurious power may arise as a result of the finite extent of 
the time series. These regions are referred to as the cone 
of influence (COI), having a temporal extent equal to the 
e-folding time of the wavelet function. For the case of the 
Morlet wavelet it has the form, = V2s = \/2P/1.03. 
This time scale is effectively the response of the wavelet 
function to noise spikes and is used in our detection crite- 
ria by requiring that oscillations have a duration greater 
than td outside the COI. This imposes a maximum period 
above which we can not accept detections, since we are 
dealing with finite time series of duration D to t ■ At this 
maximum period it follows that P max = 1.03D to t/^V^ 
— for set 1 D tot is 2032 s and P ma x is 493 s, while D to t 
is 7060 s and P ma x is 1714 s for set 2. 

The form of automated wavelet analysis carried 
out on these data has previously been presented in 
iMcAteer et al.l()2004|) . although some revisions have been 
carried out for our work. Firstly, no filtering was applied 
to data set 1 because the duration of observation pre- 
sented here (34 min) is small compared to that of the 
TRACE spacecraft orbital period (^90 min). Therefore, 
long-period variations due to changes in Earth upper- 
atmospheric transmission are of small amplitude when 
compared to shorter periods of solar origin. Secondly, 
the more stringent statistical criterion of oscillations be- 
ing at least 99% confident is employed rather than pre- 
vious values of 95%. Finally, all power in the transform 
is retained, even that in the COI. Any detection having 
power above the 99% level for at least outside the COI 
has the entire duration recorded, including any portion 
which falls in the COI. This modification helps somewhat 
to circumvent previous limits on durations of detections. 
Not only does the upper limit to the possible number of 
cycles decrease with increasing oscillatory period (since 
analyzing finite duration time series), but removing the 
COI decreases this upper limit more drastically as the 
COI both increases with period and occurs at the start 
and end of the wavelet transform. After a complete anal- 
ysis, carried out on the time series of each spatial pixel 
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Fig. 3. — Fourier phase differences (ordinate) between 2420 
pairs of TRACE 1700 A and 1600 A time series as a func- 
tion of oscillatory frequency (abscissa) . Larger symbols and 
darker shading depicts increased weighting significance (i.e., 
increased Fourier cross-power and coherence) . The phase dif- 
ference gradient (M^\ solid line) was calculated by weighted 
least-squares fitting of a first-order polynomial to all points 
in the range 2.0 — 5.4 mHz, indicated by vertical dotted lines. 

in all three TRACE passbands, the output comprises the 
duration in oscillatory cycles and the period of every de- 
tection. 

3.2. Fourier Phase Differences 

A Fourier phase difference analysis was used to fur- 
ther study possible propagation characteristics of waves 
between the co-spatial 1700 A and 1600 A time se- 
ries. These passbands were chosen b ased on the proxim- 
ity of their peak formation heights (Foss um fc Carlssonl 
2005b). This form of Fourier analysis is a technique 
which allows investigation of the difference in waveforms 
between two time series. Using the phase information 
contained within the complex Fourier transform, vi a the 
equations discussed in detail by Krii ger et alJ (|2Q0 lf> . the 
difference in cyclic phase, A<f>, can be determined for 
each Fourier frequency component, v, with the quality 
of the values represented by the phase coherence. The 
pairs of 1700 A and 1600 A time series in this study are 
extracted from the same (x, y) pixel location, but the 
signals remain separated in the direction normal to the 
solar surface. In this scenario, phase differences can be 
interpreted as delays caused by the finite propagation 
speed of waves traveling between the formation heights. 
The 1700 A data were linearly interpolated to the 1600 A 
(equally-separated) times to compensate for timing off- 
sets that result from the TRACE UV data acquisition 
(§HJ> ■ These offsets require removal as they produce drifts 
in ph ase difference spectra with frequency ijKriiger et al.l 

At those frequencies over which a wave has a real solu- 
tion to its dispersion relation, and may thus propagate, 
phase difference spectra between emission signals formed 
at two separate heights are expected to have the form, 



Ah 
— i 

Vph 



(1) 



where Ah is the separation distance of the emitting 
regions, v p h is the phase speed of the wave, and 
v is the oscillatio n frequ ency. In a method similar 
to IMcIntosh et al.1 l)2003j) . phase difference gradients 



(Ma^ ~ Ah/ Vph) were extracted over the frequency 
range 2.0 — 5.4 mHz (periods 185 — 500 s). This range 
was chosen because it encompasses the majority of our 
oscillation detections (§ 14.1(1 . As shown in Figure^ phase 
difference gradients were achieved by applying a weighted 
least-squares fit of a first-order polynomial to the phase 
difference spectra obtained from all time series originat- 
ing from pixels in a particular range of field strengths. 
The contribution of each phase difference point to the 
fit was weighted by both the Fourier cross-power and 
coherence values. 

4. RESULTS & DISCUSSION 

Wavelet results are presented only for set 1 to avoid 
repetition, since set 2 displays similar behaviour with 
longer-period detections due to the increased duration of 
observation and more short-period detections from flare 
surges in the time series. However, Fourier results from 
both data sets are provided and discussed. 

4.1. Wavelet Period Distributions 

Detections in each mdi field strength interval were sam- 
pled in terms of their period with uniform 20 s binning. 
The resulting 1700 A, 1600 A, and 1550 A PDFs are 
shown in Figure^] where the individual distributions are 
self-normalized to their peak number of detections to 
highlight changes between PDFs by compensating for the 
decreased number statistics at increased field strengths 
(Figure 0) ■ Fitting profiles which consisted of a single 
Gaussian and a linearly increasing first-order polyno- 
mial (as most distributions showed a tail of detections 
at high periods) were applied to each PDF. Although 
the distributions are not ideally Gaussian, the variation 
of fitted centroid period (P c ; displayed as central bar 
marks in Figure0|and detailed in Table[2J with mdi field 
strength and UV passband formation height should at 
least qualitatively describe these data. Bars extend over 
the ±lcr width of the Gaussian fit in Figure^] represent- 
ing PDF spread and not error in the fit. For each UV 
passband the symbols are plotted at the mean absolute 
field strength of the pixels which show oscillations in that 
passband. For a particular range of mdi field strengths 
the mean value will vary between each of the three UV 
passbands, a direct result of the differing pixel numbers 
which contribute to the different UV passband oscillation 
distributions (lower panel of Figure 

Period PDFs in the range |25 — 100| G shift to higher 
values of P c when moving to larger MDI field strengths. 
This progression of periods in the low chr omosphere from 
4- to 5- min agrees with the findings of lLawrence et ail 
( 2003) who analyzed 3 A-wide passband observations of 
Ca II K line intensity in the quiet-Sun (inter)network. 
Although the K2 portion of the Ca II K line pro- 
file is formed at least somewhat above the tempera- 
ture minimum (i.e., above these TRACE data), close 
agreement may still be expected with this work be- 
cause the Ki wing emission also contained in their broad 
passband is formed a round the temperature minimum 
l(Vernazza et al.l fl981^l . Lawrence et al. attribute the 
period change to their passband being formed below 
the f3 = 1 level at low field strengths and above it for 
higher fields, without offering a physical reason for the 
change. We draw attention to the fact that the periods 
detected are incompatible with standard MHD theory, 



TRACE UV Wave Dependence on Magnetic Field 



5 




TABLE 2 

Magnetic Field Strength Intervals and Fit Parameters to 
Wavelet Period Distributions of 2004 July 03 (Set 1) 
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Fig. 4.— Frequency distributions of TRACE 1700 A (lower), 
1600 A (middle), and 1550 A (upper) wavelet periods (ordi- 
nate) as a function of the mdi photospheric magnetic flux 
magnitude (abscissa) . Distributions are normalized to their 
peak number of detections. Symbols mark centroid positions 
of Gaussian fits to each period distribution, with bars ex- 
tending over the ±1<t width (detailed in Tabled- Points are 
overlaid at the mean mdi value of pixels showing detections 
in that field strength range. The dashed line marks the non- 
magnetic value of P ac for 9 — at the temperature minimum. 



where waves can only propagate at periods below the 
acoustic cutoff (P ac oc VT; T is temperature). As such, 
waves up to ~5 min period can propagate into the photo- 
sphere (e.g., p modes) but the temperature decrease with 
height in the low chromosphere acts to decrease the cut- 
off period and previously propagating long-period waves 
become evanescent. 

We extend the hypothesis of lLawrence et al.1 l)2003ft by 
proposing that wave propagation may be permitted and 
the period behaviour explained using both the direct de- 
pendence that P ar has on the magnetic field strength 
(iBel fc Meinl and the inclination, 9, which arises 

through the effective grav ity, q (P ar oc 1/q = 1/qn cos 9, 
where g = 274 m s" 1 : iDe Pontieu et alJ 12004). As 
shown by lLawrence et alJ l)2003ft . regions of increased 
magnetic field experience the (3 = 1 level at lower al- 
titudes. This intuitively suggests that the "magnetic 
canopy" , where the field spreads nearly horizontally, will 
be formed lower. If for low field strengths a HOF lies 
significantly below the "canopy" level it will experience 
near-vertical field lines. However, as the photospheric 
flux increases, the (3 = 1 and "canopy" levels will be 
pulled down toward the HOF such that it samples in- 
creasingly inclined fields. The larger inclination acts 
to reduce the effective gravity, which in turn increases 
the cutoff peri od allowing for "tunne ling" of evanescent 
wave modes l|De Pontieu et all 12004ft . Over the range 
1 22 — 75 1 G, increasing values of P c are recorded when 
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moving upward through the UV passband HOFs, agree- 
ing with this picture of sampling field lines of increasing 
inclination with height. 

All three UV passbands subsequently show decreasing 
centroid periods when moving into regions of increased 
magneti c field (|125 — 5 5 0| G) . This differs from the re- 
sults of lLawrence et all l|2003ft where the peak wavelet 
power was observed at a fairly constant 5-min period for 
field strengths in the range 1 200 — 500 1 G, although a 
marginal downturn in period may exist at their highest 
fields. However, it sh ould be noted that th e solar en- 
vironment studied by lLawrence et alJ (J2003) contained 
neither plage nor sunspot regions which contribute to 
those field strengths in this work. A more appropriate 
comparis on can be made with data set 2 previously pre- 
sented in .Mcintosh & Smili^ l)2004ft which is our set 2. 
Although not discussed, comparing the detected wavelet 
frequency in their Fig. 5 with the magnetic field shows 
the same variation: short period oscillations in the inter- 
network, increasing periods in areas of network and outer 
plage, and a decrease in period in those plage regions cor- 
responding to the strongest plage magnetic fields. 

The mdi intervals containing all but the largest field 
strengths (|650 - 950| G) exhibit double-peaked PDFs. 
This makes qualitative sense as a sunspot lies in the FOV. 
As such, the lower-perio d components co ntain the ~3- 
min umbral oscillations ijGiovanellil 1 1 9721 more distinct 
in data set 2 for fields > 1 1200 1 G where periods peak at 
160— 180 s), while components at hig her periods contain 
the ~5-min penumbral os cillations ijBeckers fc Schult3 
Il972t iMoore fc Tanal 11975ft . It should be noted that 
period PDFs for all passbands do not appear double 
peaked at the highest fields (> |950| G), despite show- 
ing decreased number statistics. The disappearance of a 
higher-period component at these largest field strengths 
supports their interpretation as penumbral oscillations 
since these fields should sample only umbral regions. 
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Fig. 5.— As Fig. E] but for TRACE 1700 A {lower), 1600 A 
(middle), and 1550 A (upper) durations of wavelet detection 
(ordinate). Symbols mark the most frequently detected du- 
ration, plotted at the mean field strength of pixels showing 
oscillations in that field strength range. 



4.2. Wavelet Duration Distributions 

An alternate way by which to characterize oscillations 
is their duration of observation, since this quantity can be 
interpreted as a measure of the energy associated with a 
wave. As such, detections were sampled in terms of their 
duration with uniform half-cycle binning. The resulting 
PDFs are shown separately in Figure El for the 1700 A, 
1600 A, and 1550 A data, where every PDF is again 
normalized to its peak number of detections. In these 
plots the most frequently detected duration is marked by 
a '+' symbol, and differing dependencies on the magnetic 
field are clearly seen. The 1700 A data have oscillations 
of duration 2.5 — 3 cycles most frequently observed across 
the majority of magnetic fields, while the 1600 A data 
also peak at 2.5 — 3 cycles over most field strengths. The 
large degree of similarity between the behaviours of these 
two sets of distributions strongly implies that both UV 
passbands are formed in essentially the same physical 
environment. However, the 1550 A durations peak at 
1.5 — 2 cycles over low field strengths ( 1 25 — 100 1 G) before 
changing to 2.5 — 3 cycles for moderate magnetic fields 
( 1 100 — 450| G). The oscillation durations then shorten 
to 2 — 2.5 cycles for nearly all higher field strengths. 

It appears that, if these oscillations are indeed waves, 
below |100| G some portion of the wave energy budget 
has been lost between the formation heights of 1600 A 
and 1550 A, since the 1550 A data exhibit durations 
shorter than both the 1700 A and 1600 A data. By con- 
trast, the wavelet detection durations are conserved in 
the stronger-field range 1 100 — 450| G which, in turn, 
indicates that the effective "boundary" previously ly- 
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Fig. 6. — Fourier phase difference gradients (Ma^; ordi- 
nate) between the 1700 A and 1600 A intensity time series 
as a function of mdi magnetic flux magnitude (abscissa) — 
triangles correspond to set 1 and squares to set 2. Symbols 
are plotted at the mean absolute field strength value, where 
horizontal bars denote the extent of the field strength range 
and vertical bars mark the ±lcr uncertainty to the weighted 
least-squares linear fits. The values obtained over low-field 
values are displayed within the inset to highlight fine detail. 

ing between the 1600 A and 1550 A HOFs has been 
lowered in altitude to below the 1700 A HOF. Incor- 
porati ng a potential-field extrapolation, lLawrence et alJ 
( 2003) interpret their change from 4- to 5- min periods 
at field strengths of ~|150| G as resulting from the low- 
ering of the /3-unity surface below the formation height 
of their Ca II K line emission. This is very close to 
our wavelet duration transition field strength of |100| G, 
while the disparity may be explained by the difference 
in passband formation heights and/or the assumption of 
the static, spatially-a veraged VAL C model atmosphere 
l)Vernazza et al.iri981|) in the extrapolation. 

4.3. Fourier Phase Speeds 

Although highly versatile and a powerful tool for 
studying oscillatory behaviour, the wavelet analysis per- 
formed up to now has not been able to distinguish 
whether the detected oscillations are the signatures of 
propagating waves in the solar atmosphere. It is for this 
reason that Fourier phase difference analysis was applied 
to the same data. The results of the analysis, outlined in 
§ 13.21 are presented in Figure El for both data sets as the 
weighted-fit values of phase difference gradient between 
the 1700 A and 1600 A time series pairings. 

If a spatially-invariant atmosphere is assumed, the ver- 
tical separation of UV HOFs will not vary when mov- 
ing into regions of increased magnetic field. As a re- 
sult, Ah is constant and Equation ^ informs us that a 
change in Fourier phase difference gradient denotes the 
inverse of the change that has occurred in the wave phase 
speed. In the past, such forms of atmosphere have been 
used to implement potential-field extrapol ations (e.g., 
IMcTntosh et ahlEfiM IT^wTence et al .1120031. Unfortu- 
nately, these act to oversimplify the complexity of the 
magnetic atmosphere and the application of 1-D models 
that take account of the differences between regions of 
differing magnetic field will yield more appropriate re- 
sults. 

From consideration of Equation ^ it is possible that 
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the initial decrease of gradient values in Figure could 
result from waves having the same phase speed but the 
HOF separations are reduced with increasing magnetic 
field. This suggests that Ma^ should continue decreas- 
ing for even stronger fields, but we detect a switch in 
behaviour to increasing Ma,/, in the range 1 125 — 550 1 G 
which requires increasing separations between the HOFs. 
This change in separation behaviour does not make qual- 
itative sense, so we proceed by assuming a spatially- 
invariant atmosphere that yields increasing wave phase 
speeds up to 1 100 — 125| G and a subsequent decrease 
in speeds up to |550| G. A decrease in TRACE UV 
phase difference gradient with in creasing magnet i c field 
has been previously observed by [Mcintosh ct al. (2003) 
with no turnaround reported. The deviation of our re- 
sults can be explained by their quiet-Sun (inter)network 
data having fields <|50| G so the (3 = 1 level was not low- 
ered below the UV HOFs, while wavelet durations imply 
that this occurs for fields >|100| G in this study. 

The observed wave speed behaviour over low- to 
moderate- magnetic fields may be interpreted as result- 
ing from a change in physical environment within our 
pixel sampling rather than merely a variation of mag- 
netic field. In regions spatially removed from the large 
flux concentrations of sunspots and plage, photosphcric 
magnetic fields are believed to be unresolved kilo-Gauss 
elements, covering a small fraction of the pixel area 
l)Litesl2002tlDominguez Cerdeha fc Kneerl2003|) . In such 
cases the magnetic filling factor, /, is much less than 
unity and pixels will contain a contribution from both 
the ambient high-/3 plasma and low— /3 flux concentra- 
tions. In the high-/? material, magnetoacoustic slow 
modes have negligible density variations and are unob- 
servable to our broad UV passbands, while the compress- 
ible fast modes will propagate between eg and ycf + v\ 
(~ cs since |3 > 1). In the flux concentrations, slow- 
magnetoacoustic waves will propagate at cs (since (3 < 1) 
while fast modes propagate between va and ycf + v 2 A . 
If these flux concentrations satisfy certain requ irements 
they may be thought of as "thin flux tubes" l)Robertsl 
1985), in which slow modes propagate at the tube speed. 
Values of 0.2 — 0.5 are typicall y used for plasma-/ 3 in the 
modeling of flux tubes fe.g.. iHasan et al.ll2003l [2005), 
yielding tube speeds of ^0.9cs and fast-mode speeds 
between ~(1.5 — 2.4)cs and ~(1.8 — 2. 6)c£. In this pic- 
ture, the lowest field strengths arise from small values 
of / so the measured speed is dominated by the "flux- 
external" sound speed, while for higher field strengths 
the filling factor increases and the contribution from the 
"flux-associated" speeds becomes more important. 

Eventually, a filling factor of unity is reached and the 
passbands are formed in an homogeneously distributed 
magnetic environment for greater field strengths. The 
change from inhomogcncous to homogeneous environ- 
ments ends the "flux tube" supported wave scenario. 
When / = 1, the magnetic field has expanded to fill 
the entire pixel area, creating a "magnetic canopy" over 
the previously "flux-external" plasma. This marks the 
transition from mixed high- and low- (3 to entirely low— /3 
environment (i.e., the (3=1 level) and is signaled by the 
change in trend of phase difference gradient at fields of 
1 100 — 125 1 G. Although kilo-Gauss magnetic elements 
have been mentioned and field strengths of 1 100 — 125 1 G 



are reported for filling factors of unity, field values are 
measured at the photosphere while / corresponds to the 
region of the chromosphere where waves are detected and 
field lines have expanded into a greater spatial area. 

Inside the subsequent homogeneous low-/? atmosphere, 
the fast-mode propagates at Va > cs, tending to 3>cs 
when (3 -C 1. The initial fast-mode speed in the ho- 
mogeneous situation (i.e., when (3 = 1) is expected to 
be at most ~\/2cs, yielding a phase difference gradi- 
ent which differs from the field-free case by a factor of 
~0.7 (within the error bars of the ratio between our 
values recorded at 1 100 - 125| G and |0 - 20| G). The 
slow mode is again quasi-longitudinal and propagates at 
cs- We note that the sound speed should not change 
between field strengths since it depends only on tem- 
perature ((xVT) and the wave behaviours reported here 
occur between the same two UV passbands, which are 
formed over the same range of temperatures irrespective 
of magnetic field. The gradual decrease of wave speeds 
in the range 1 125 — 550 1 G is taken as a signature of 
both modes existing within our signal, with the mea- 
sured value of Ma^ resulting from averaging between 
the two speeds. The increase of phase difference gradi- 
ent with increasing magnetic field towards a similar value 
as detected over field strengths in the small filling factor 
range of |0 — 20 1 G implies that the wave speeds are 
tending toward the sound speed, and slow-mode waves 
become dominant in our stronger field regions of plage 
and umbrae. An implication of slow-mode dominance 
in regions of stronger field is that slow mode s propagate 
highe r into the atmosphere than fast modes l)Osterbrockl 
I1961D — the "magnetic canopy" occurs at lower altitudes 
when moving to higher field strengths, the UV passbands 
are formed progressively further above the /3-unity level, 
and we detect only the speeds of modes which propa- 
gate higher into the atmosphere. The disclosure of these 
waves as slow modes p rovides a direct, and ne cessary, 
link to coronal studies l)De Moortel et all 120021) . where 
propagating intensity variations at 3- and 5- min periods 
are interpreted as slow-magnetoacoustic waves in loops 
anchored in umbrae and plage, respectively. We note a 
disparity between the abundance of wave signatures in 
the low chromosphere and the sparse detection of slow- 
mode associated pheno mena in the transition region 
llDe Pontieu et all 120041) and corona i|De Moortel et alJ 
2002). Some mechanism is clearly required to quench 
most of these lower atmospheric waves with one possibil- 
ity being slow-m ode driven Pedersen current dissipation 
( Goodmanl l2000j) . theorised to be highly efficient in the 
chromosphere. The non-uniformity observed in the up- 
per atmosphere may highlight other selection effects at 
work on wave transmission, possibly based on p-mode in- 
terference with complex field orientations. Field inclina- 
tions affected by p-modes may also help explain the tem- 
poral intermittency of waves by only providing suitable 
propagation conditions over short, repeating timescales. 

It should be noted that some penumbral pixels may 
be misinterpreted as weaker field regions from of the 
LOS nature of mdi and the existence of highly-inclined, 
outer-penumbral magnetic fields, e.g., fiel d strengths of 
1800 - 1000 1 G inclined 80° to the vertical ijSolanki et alJ 
1992) will be recorded as 1 140 — 175 1 G. However, pixels 
misclassified as such will not greatly modify the phase 
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difference gradients, as the true weak-to-moderate field 
pixels (| 25 — 550 1 G) dominate the number statistics. In- 
stead they contribute to greater error-bar magnitudes 
at higher field strengths, alongside the decrease in pixel 
numbers. Larger errors associated with gradient val- 
ues above |550| G may arise from averaging between 
the sunspot ~3- and 5- min distributions, if propagat- 
ing at different speeds. Decreased gradient values over 
1 550 — 950 1 G from those at fields both above and below 
complements this interpretation, as period PDFs exhibit 
widening toward and the detection of 5-min penumbral 
oscillations over this range. Hence, by causing decreased 
phase difference gradients, penumbral waves appear to 
either travel faster than those in umbrae and plage or ex- 
perience a reduction in the HOF separations, distinctly 
possible for a penumbral atmosp here stratified along f ield 
lines at inclinations of 40 - 70° ijMathew et al.ll2003jh 

5. CONCLUSIONS 

Chromospheric intensity variations have been studied 
using both time-localized wavelet and Fourier analysis. 
Investigation of the oscillation behaviours has identified 
three regimes with distinctly differing properties. 

The increase of wavelet period with MDI field over 
1 25 — 100 1 G marks the modification of the magnetoacous- 
tic cutoff period by field inclination effects. This arises 
from the "magnetic canopy" being pulled progressively 
closer to the sampled HOFs so that they experience more 
inclined fields. Wavelet durations show a decrease in the 
1550 A data over those in the 1600 A and 1700 A data, 
indicating a loss of wave energy. This can be explained 
by waves undergoi ng mode conversion or reflection at the 
/3-unity surface ijBogdan et alJl2003l) . Fourier analysis 
indicates that wave phase speeds increase with magnetic 
field between the 1700 A and 1600 A formation heights. 
It is proposed here that this is due to an increase in mag- 
netic filling factor within the pixel sampling, yielding a 
greater contribution to the observed wave speed from the 
low-/? "flux-associated" fast-mode speed over the slower 
sound speed of the surrounding high-/? plasma. 

The similarity of wavelet periods and durations over 
field strengths in the range 1 100 — 450 1 G implies that 
waves in this regime do not experience a drastic change 
in their environmental topology as evidenced by the 
weaker-field case. The /3-unity surface is believed to lie 
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below the formation height of 1700 A emission, no tran- 
sitional interface exists between the passbands, and wave 
energy (i.e., wavelet duration) is conserved with altitude. 
In contrast to the weak-field regime wave speeds de- 
crease, which may be a signature of a change from dom- 
inant fast-mode magnetoacoustic waves to slow-mode 
waves when moving from network into plage. 

Wavelet period distributions broaden and display a 
doubly-distributed nature over 1 550 — 950 1 G, marking 
the sampling of sunspot regions. Interpretation of the 
5-min component as being of penumbral origin is given 
strong support by the virtual disappearance of these 
higher-period components when sampling field strengths 
>|950| G (i.e., umbrae). The decrease of phase differ- 
ence gradient values over the mixed penumbral/umbral 
field strength range (|550 — 950| G) suggests that either 
penumbral waves travel at greater speeds than those in 
plage and umbral regions or the formation heights of the 
1700 A and 1600 A passbands are closer in altitude (e.g., 
because of inclined fields) in the penumbra. 

In the future, we will be able to better determine 
the intricate relationship between magnetic topology and 
propagating waves using high-resolution, vector mag- 
netic field information. Such data will soon be obtained 
by the full-Stokes spectropolarimctcr in the instrument 
suite of the Solar Optical Telescope (sot) Focal Plane 
Package (fpp) on board Solar-B. Coordinated observa- 
tions between this instrument, TRACE, and the multi- 
wavelength tunable filters also contained in the SOT / fpp 
will yield diffraction-limited, high-cadence imaging at 
multiple heights through the outer Solar atmosphere un- 
paralleled by previous observations. The combination of 
these forms of data with a range of atmospheric models 
taking into account the differences in atmospheric strati- 
fication for differing magnetic fields will yield a powerful 
diagnostic suite for studies of solar wave phenomena. 
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